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ABSTRACT 

The optimal dumping of organic waste and the artifical aeration 
of a river in order to maintain adequate levels of dissolved oxygen in 
the water are considered in this thesis. The problem is formulated as 
an optimal control problem with several performance indices. Available 
results on ren ene optimization theory are used in the present study. 
Optimum profiles of BOD (Biochemical oxygen demand) and DO (dissolved 
oxygen) are obtained for several objective functionals. These functionals 
are expressions relating the chosen performance indices. In previous work 
in this area, the approach has been to use a single performance index. 
However, since a river has to fulfil many requirements,some often 
contradictory, the multi-cost approach used in this thesis has a decided 


advantage. 
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CHAPTER (1) 
INTRODUCTION 
1.1 Man and his Environment 

It is becoming increasingly evident that the environment is not 
limitless.* In some areas, pollutants and associated destructive agents 
are affecting the health, recreational life and the opportunity to work 
of each and every man. Air is a continuous requirement of man while 
water is consumed selectively in different forms along with food from 
the soil during certain times of the day. These basic requirements for 
life are endangered, since the ecosystem is being altered by man's 
activities.** Thus engineers are becoming concerned about this problem 
and are attempting to do their share to correct the wrongs which have 
resulted from man's previous indifference to the ecosystem. 

The soil has been relegated to the role of the recipient of 
many products resulting from man's activities such as solid refuse disposal 
from both urban and rural sources, chemical agents used in the production 
of food, and contamination froin particles settling out of the atmosphere. 
It has been recently discovered that the use of pesticides has caused 
harmful residues to be concentrated in the food chain. The pollution of 


air is caused by burning of fuels for transportation, heat and power. 


* Environment: a continually changing complex of all conditions and 
influences interacting with an organism. 


*k Ecosystem: a relationship between living things and their nonliving 
but supporting environment. 
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Sources of air pollution emit some or all of the following: sulfur dioxide, 
particle matter, carbon monoxide, photochemical oxidants, hydrocarbons 
and nitrogen dioxide [2].* It has been speculated that the oxygen - 
carbon dioxide cycle in nature might be upset due to these pollutants. 
In the case of water, municipalities and industries use lakes, rivers 
and oceans for the dilution of many liquid wastes which they generate. 
By using nearby rivers or lakes for the disposal of liquid wastes 
industries such as pulp and paper have introduced toxic substances like 
mercury into the water. This has resulted in the recall of some food 
products and the banning of others, for example certain types of fish, 
from the market place. 
1.2 Water Pollution 

The general quality of water in lakes, rivers and oceans, 
depends mainly on mereacatedica’ factors, municipal uses and industrial 
uses in their vicinity. The deposition of certain dissolved gases and 
particles from the atmosphere and the washing of organic matter, suspended 
solids, and chemical residues from the surface of the earth are the main 
ingredients of meteorological poliution of surface waters. Municipal 
discharges contain organic matter, infectious agents, suspended solids 
and detergents, while industrial discharges are characterized by chemical 


residues, metal ions, heat, organic and inorganic matters. 


* Numbers in rectangular brackets refer to references listed under 
Bibliography at the end of this thesis. 
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We shall discuss one aspect of this problem, namely the addition 
of organic waste to a river. However before this physical problem is 
defined, it seems advantageous to describe what is known as the “organic 
cycle of nature" [26]. 

Organic Cycle 

When food is consumed, the outcome is the liberation of energy 
and the replacement of worn out cells. Nitrogen, phosphorous and other 
elements are used to repair the deteriorated cells. The replaced 
materials are discharged through the human system as waste products. The 
energy content of these waste products are lowered further by bacteria 
acting on it. Each genus of bacteria in turn lives on the residues or 
sewage of the higher order species, until such time, no other organism can 
derive any further benefit by this process. Biochemical degradation refers 
to the above reaction brought about by biological agents to lower the 
energy contained in the molecular structure or to simplify the molecular 
structure. The stabilized end product is the fertilizer from which crops 
derive paca to complete the organic cycle. 

The above organic growth and decay occur under two different 
sets of conditions. One proceeds when a sufficient quantity of free 
atmospheric oxygen is available for biodegradation by organisms in the 
water. This process is referred to as “aerobic decomposition". The 
other process is "anaerobic decomposition" in which organisms in the 
water release oxygen from other chemical compounds for the stabilization 
of organic matter. The energy released by the bacteria under the aerobic 


conditions is about thirty times [26] that is available to the bacteria 
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under anaerobic condition. This accounts for faster stabilization by the 
former rather than by the latter. In addition, the anaerobic end products 
are not stable. Thus they are still subject to subsequent aerobic 
digestion when the conditions are right. For the complete description 

of the two processes refer to figure 1 [26]. 

From the municipal sewage the return water, carrying man's 
ground garbage, human feces, and urine, will be subject to the right hand 
half of the different cycles when the unstable organic compounds are 
rendered into the river. The amounts of these unstable compounds can be 
measured by a quantity known as the biochemical oxygen demand (BOD) which 
is usually defined as the amount of oxygen required by bacteria while 
Stabilizing decomposable organic matter under aerobic conditions. The 
ultimate BOD of the added sewage is the total amount of oxygen required 
to reduce the organic matter to stable compounds. Sometimes the oxygen 
demand proceeds at rates slower than natural reaeration of the receiving 
water, while at other times it far exceeds these rates. These levels of 
the dissolved oxygen (DO) are derived directly from the air or indirectly 
through photosynthetic processes of aquatic plants. It should be obvious 
that BOD and DO are the prime indicators of the quality of water in a 
river after organic sewage is dumped into it. 

For perspective the BOD from all municipal sources throughout 
Canada for the twentieth century is presented in figure 2 [12]. It is 
an increasing curve since the amount depends on the population and the 
standard of living of the country. The present state of technology 


can limit the final concentration of BOD added to a river by secondary 
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FIGURE 1: Nitrogen, Carbon, and Sulfur Cycle in Aerobic and 
Anaerobic Decomposition 
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sewage treatment which achieves a ninety percent removal rate and by 
tertiary sewage treatment which removes an additional five to ten percent. 
Even with the tertiary system the concentration of the BOD dumped into a 
river is one to four and half percent of the initial value received by a 
treatment plant. 
1.3 Application of Modern Control Theory to Water Pollution 

In the past few years there has been a number of papers published 
on the application of control/system theory to the problem of adding waste 
to a river [1,6,7,8,9,22,23,24]. This research has been concerned with 
the relationship between the concentration of DO and BOD caused by 
decomposing organic matter in a waterway. The problem is formulated as 
an optimal control problem in which the DO and BOD are used as controls 
to minimize a performance criterion according to calculus of variation 
techniques. In several of the papers [1,6,22,23] the objective functional 
is formulated either as a single index or as a linear combination of several 
performance indices using weighting factors to relate them. These weight- 
ing factors are arbitrarily specified a priori in the papers by Tarassov, 
Perlis and Davidson [23] and by Fan, Nadkarni and Erickson [6]. This 
approach has a drawback since it implies that a known relationship exists 
between the different performance indices used by them, which may not 


always be the case. 


1.4 Scope of Thesis 


The intent of research reported in this thesis is to determine 
the optimum dumping control, i.e. the addition of organic waste, and 


the optimum aeration control, i.e. the addition of oxygen by artificial 
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FIGURE 2: BOD from Municipal Sources throughout Canada 


or 


means, in terms of the weighting factors for all possible objective 
functionals relating the selected performance indices. In this manner 
more insight can be gained into the behaviour of the system. 

In chapter two, a lumped parameter model in state space form 
is derived. The concentration of BOD and the deviation of the concentration 
of the actual DO from its saturation value are chosen as the state 
variables. This model is compared with a bilinear model which is also 
derived in this chapter. 

Chapter three discusses the criteria for optimum water quality. 
Performance criteria and objective functionals are formulated. 

In chapter four, optimization of water quality is carried out 
using the multi-cosSt optimization approach developed by Salama and 
Gourishankar [20]. 

In chapter five, the results are summarized and conclusions 


stated. 
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CHAPTER (2) 
MATHEMATICAL MODEL 

2.1 Physical Problem 

The system to be modeled consists of a section of a river with 
a continuously-treated municipal waste discharged into it. (See figure 3). 
The addition of artificial aeration along the section can be included, if 
necessary. Starting with a distributed parameter model used by earlier 
investigators, it will be shown that by making appropriate simplifying 
assumptions, a lumped parameter linear model in state space form consisting 
of two first order equations will be derived. A bilinear model will 
also be described. 
2.2 Distributed Parameter Model 

From the previous chapter we have determined that it would be 
advantageous if we could determine the demand and the supply of DO at 
any point in a section of a river. In general, the movement and reactions 
of pollutants being transported in a river, can be described by the 
conventional diffusion equation [19] having the form 

Zz Z 526 
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y 
where c is the concentration of the pollutant 
t is the time 
x is the longitudinal distance along the river 
y is the width of the river 
z is the depth of the river 


D. is the dispersion coefficient of the river (j=x,y,z) 


r(c) is the rate of pollutant decay 
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Municipal Waste 
Treatment Plan: 


Effluent Outler 


Continuous Discharge of BOD 
within the River 


el Pa Osan Nes cnccen asa Ye — 
——— Or ee me ti) pe ye 
TN 


FIGURE 


3 


Schematic of River Section 
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V(x,t) is the velocity of the river. 

If we recognize that the width and depth are small compared to 
the length under consideration, then it can be assumed complete mixing 
is achieved within a short distance downstream. Now the terms involving 
y and z disappear and the concentration of the mixture depends on the 
downstream distance from the beginning of the river section. The effects 
of dispersion along the length of the river are small compared to the 
rate of transportation downstream due to the velocity of the river. 


Equation (2-1) becomes [5,11,23], 


dee = OCS - 
ahs -V (x,t) = ECe) (2-2) 


Using the above general equation, two equations can be 
written to describe the level of poilution due to organic waste. The 
first is concerned with the BOD due to the oxygen required by the bacteria 
to stabilize the decomposable organic waste under aerobic conditions. The 
second equation gives the DO derived from the air and the aquatic plants. 


The equations are 


OB(x,t) 


at = “K, (x,t) - V(x,t) oe 9 (2-3) 
OD( x,t) : aD(x,t) 
Ot = “K,(,t) B(x,t) V Gx, t) nae + 


KGs,t) [D, (x,t) - D&x,t)] (2-4) 
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where D is the DO concentration 
B is the BOD concentration 
Ko is the reaeration coefficient 
K. is the BOD removal coefficient 
Ka is the deoxygenation coefficient 


De is the saturation level of the DO. 


Next we consider a section of a river which has uniform 
properties. In other words, the above coefficients are constant with 


respect to time and distance. Equations (2-3) and (2-4) become 


SBCx,f) = -K. Bx,t) - V(x) oBCx,t) (2-5) 
SPs) = -K, D(x,t) - Kae B(x,t) + a Ne 
V(x) oe) (2-6) 


Equations (2-5) and (2-6) constitute a distributed parameter model of 
the waste disposal system. This model has been used by some of the 
earlier investigators [18,22,23]. 
2.3 Lumped Parameter Model 

If the concentration of BOD from the dumping source is constant 
with respect to time, we get the following lumped parameter model of the 


system [1,4,9,18,20], 


aw (0) (2-7) 
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describing the concentration of the BOD and DO for a section of a river. 
As can be seen from equations (2-7) and (2-8), we have not 
taken into consideration the addition of organic wastes in terms of BOD. 
Since the levels of BOD and DO along a section depend on the inflow of 
organic waste, we define u, (x) to be the control function representing 


this input. Now the mathematical model becomes 


K u, (x) 
BOO oa Ae 1 — 
age. Orta Gath xhotesy (2-9) 
K K 
aD (x) eed ies > = 
Se So SG HEN (2-10) 


Let us now cast the model in matrix form using the following 


state variables representation, 
x, Cy) = B(x) 
x, (y) = D. - D(x) 


where y replaces x, the distance downstream for a section of 


a river, to avoid confusion in the notation. 


Thus, 
ky) sty 0} (x) | 
= | + | | [u, Gy] (2=11) 
x (J Le Ail) Lod 


Initial Conditions: x, (0) =x and x» (0) = x 


10 20 
where A = K,/V, B = K/V, C= K4/V and u, Cy) = u(x)/V. 
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A block diagram representation of (2-11) is shown in figure 4. 
2.4 Bilinear Model 

A bilinear model is now set up which takes into account the 
effect of DO on the rate at which BOD declines along a section of a 


river. The form of the bilinear model is 


Bx) --= D(x) + B(x) (2-12) 
OG) ~ -B aq) + p, - 0@) (2-13) 


The only difference between the two models lies in the equation 
for the rate of change of BOD. The constant Ke the BOD removal coeffi- 
cient, has now been replaced by the expression K D(x), a constant 
multiplied by the level of DO along a river. This system could be 
representative of the aerobic decomposition since the BOD level would 
remain constant when there is no free oxygen present in a section of a 
river. Then another relationship would have to be developed for the 
case of anaerobic decomposition. 

Both systems, the lumped parameter model and the bilinear model, 
were solved with the aid of an analog computer (See Appendix I) under 
identical conditions, while the values of K.. and K were varied over a 
specified range. The results are shown in figures 5 and 6. The shapes 
are similar for both models which leads to the following relationship 


between the two coefficients of BOD removal, 


K 
Kowal. ne (2-14) 
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FIGURE 4: Block Diagram of System 
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when the errors due to the multiplier of the analog computer are ignored. 
Equation (2-14) is verified by the analog computer in figure 7 and by a 
ditigal computer using the Runge-Kutta Method (See Appendix ITI). 

2.5 Comparison of the Different Models 

The mathematical model used by Fan, Nadkarni, and Erickson [6] 
in their investigation is the same as the lumped parameter model derived 
here, except for the dispersion term which is disregarded. For our case 
the major mechanism of transportation for a pollutant is the velocity of 
a river rather than the much slower physical process of dispersion through 
a liquid. 

The difference between the lumped and the distributed parameter 
models is that the former allows for the determination of the optimum 
controls in terms of only distance while the latter can give the optimum 
controls in terms of time and distrance. In practice however, a control 
with only one degree of freedom is easier to implement than one with two 
degrees of freedom. 

We shall conclude this chapter with a discussion of the 
controllability of the model derived. If the system is completely state- 
controllable then for any ee each initial state x(y,) can be transferred 
to any final state x(¥,) in a finite distance Ye which is greater than 


or equal to ve Let us outline the fundamental theorem concerning this, 


Theorem I[17]: A continuous system described by 
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where x is an n-dimensional vector 


u is an r-dimensional vector 
A is an nxn matrix 
B is an nxr matrix. 


This is completely state-controllable if and only if the composite nxnr 


matrix, 


aS) (ore Gat tallle wale 


we have 


Returning to our system to check the rank of the matrix ie 


hk -K_/V 
iy 
Ege BGA’ = 
0 +K,/V 
1 a a 0 
where B = and A= 
0 +K,/V ae oy A) 


The rank of the matrix P is two, hence the system is completely 


state-controllable. It may be noted at this time, that in case the 


coefficients K. and Ka take on a range of values, the rank of matrix P 


remains constant except for Ka equal to zero. If this situation did hold 


true then the biochemical oxygen demand would have no effect on the 


dissolved oxygen level in a section of a river. This would be a 


contradiction with regards to the formulation of the problem. 


7 e 
> : , 


(a baa iy 


it A i A te 


a od rT 
Vipin 
i Woe 0 
: 4] 
a. 
W\ a ‘ 


(Letatguos ab majeye ol? added ae eb seb ei anne ee ee 
ot eas> ak 4ed9 aks sidaife 


9 wtztem to Aner any aint 


bkon bX noBcydia ekdy 1 tug 


: ee | 


eecce, t bas beorseay > nape Rh 


24%, 


CHAPTER (3) 
CRITERIA FOR OPTIMUM WATER QUALITY 
3.1 General Objectives of Water Quality Management 

In general, all physical problems have the feature that, while 
at least one solution exists, there are usually an infinite number of 
solutions possible. The aim of the optimization becomes the selection, 
out of the multiplicity of potential solutions, of that solution which 
is the best with respect to some well-defined index of performance. 
This is true of water quality management problems also. In the past a 
number of papers have defined goals to be met, in this regard Bramhall 
and Mills [3] for instance have studied the effect of public policies 
on the benefits and costs to achieve desired water quality. Brown and 
Mar [4] have determined optimum standards to be measured against alterna- 
tive water quality management policies. 

To establish the criteria for optimum water quality the 
appropriate performance indices to generate the necessary objective 
functionals have to be selected for the system under consideration. 
Many performance indices could be examined to guide the determination 
of the best conditions under which an organic waste disposal system 
should be operated. In practice however it is not an easy task to 
specify or formulate performance indices out of their technical and 
economical contexts. The reason for this is the freedom of choice 
involved and the influence this choice has on other parameters in the 
problem. One basic performance index is the cost of BOD removal. 
Liebman and Lynn /14] have used this performance index subject to DO 


constraints in their study of water quality management of the Willamette 
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river in Oregon, U.S.A. That other choices of performance indices are 


possible can be demonstrated by a brief review of some more of the 
earlier investigators in the field of water resource management. 
Fan, Nadkarni and Erickson [6] have divided a river into 


several segments and have proposed a cost benefit function on a segmented 


basis of the form: 


n 
S = ‘ie (0,Co; 4 BC. a BF) o-)) 
i=1 


where Coy is the cost of waste treatment located at the end 
of the ith segment of the river under consideration 
for water quality management 
Ca is the cost of water treatment located at the end 
of the ith segment 
BF. is the benefit function for the ith segment 
a, and B. is the nonexistence (0) or the existence (1) of 

the ith facility. 

The cost of waste treatment depends on several factors such as 
plant size, type of pollutant, degree of treatment, etc. Fan et al have 
suggested the typical curves for BOD treatment as shown in figure 8(a). 
For purposes of mathematical analysis and optimization a polynomial function 
which fits the curve is generated to give one performance index. 
Similarly, another polynomial function is fitted to the data for the 


total cost of purifying water for consumption verses the level of BOD 


in the water coming into the plant for another index (See figure 8(b)). 
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The remaining term BF, in equation (3-1) reflects the price 
society is willing to pay to maintain a certain quality of water. While 
exact values for such a function are almost impossible to estimate, 
qualitive descriptions of the benefit function are possible. Brown 
and Mar [4] have suggested typical curves which are shown in figure 8(c). 
These curves are divided into three types, "(1) those that are independent 
of water quality, such as navigation, (2) those that decrease linearly 
with water quality, such as direct water use for consumption, and (3) 
those that decrease at a given threshold of water quality, such as 
aesthetic and recreational benefits", The mathematical representation 


of the benefit function is formed as the following: 


BF, = Pp (fhe L,) + v, exp [-w, (L,/L,;)] (3-2) 


i Eee AL 
where Ps is a benefit independent of water quality 
qs is the gradient for benefits diminishing linearly 
with water quality 
L is the maximum BOD concentration in the ith segment 
L is the BOD concentration in the ith segment 
is a parameter for recreational and aesthetic benefits 
L_. is a parameter for recreational and aesthetic benefits. 
Approaching the water quality management problem from an 
economic point of view, Boyd [1] determined a water quality management 
policy which seeks to maximize the value of a river's pollutant disposal 
and water quality services with respect to limiting pollutant discharges 


and/or increasing the river's waste assimilation capacity. The allocation 


of these services depends on the maximization of the net benefit, defined 
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= i - Ge oe Oe (3-3) 
where B is the net benefit derived from the river's services 
Ze is the pollutant disposal services consumed at 
location i 
Q, is the water quality services enjoyed at location j 
K is the waste assimilation capacity investment 
C(K) is the cost of K 
X is the Lagrangian multiplier 
FE; is the marginal product of Ze evaluated at equilibrium 
pot is the marginal product of Q, evaluated at equilibrium. 
To determine the specific relationship for the above performance 
indices, Boyd used equations (2-7) and (2-8) along with experimental data 
from the Potomac estuary below Washington, D.C.. The water quality 
services enjoyed, Q, is given by the difference between the saturation 
level of DO and the DO deficit which is the DO level of the river. By 
using table 1, the waste assimilation capacity, K, for two types of 
reoxygenation devices give the cost, C(K), that would be incurred to 
offset given levels of pollutant discharge to maintain a minimum monthly 
mean of 4mg/% of dissolved oxygen. Another set of data from the Polomac 
estuary is presented in table 2, which gives the flow requirements for 


given oxygen targets, waste loads and temperatures. From this table, 


the pollutant disposal services consumed, Z, is considered to be the 
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waste load added to the river while the flow rate, F (flow augmentation 
plus natural flow), is used to maintain a certain minimum DO concentration. 
To determine the marginal products, Pe and eae the concentration of DO, 

Q, is related to the flow rate, F, of the river then using the following 


they can be derived, 


0Z/0F BZ (3-4) 


3.2 Performance Criteria for Distributed Parameter Model 

Using a distributed parameter model given by equations (2-5) 
and (2-6), Tarassov, Perlis and Davidson [23] applied optimization 
theory for determining an artificial in-stream aeration policy for 


polluted rivers. They used an objective functional of the form: 


iE x 
i f D 
J = Cy [ 4, [6 .00%= CLE xX)y dxdt 
O re) 


t x 

f£ fe f£ D ee a 

+ sf J [ettye i sdxde (3-5) 
ie) Oo 


where C(t,x) is the actual DO profile in a section of a river 
u(t,x) is the magnitude of the distributed artificial 
aeration from the in-stream aerator 


cy and c, are the weighting factors. 


The objective functional is a linear combination of two 
quadratic performance indices. The first performance index represents 


the difference between a specified value of DO level and the actual 
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DO profile. Minimization of this index yall predent the excessive 
deviation of the actual DO from the value of 6.00 mg/litre which is 
an assumed water quality standard to maintain the fish life and the 
aerobic conditions within a river. The second quadratic performance 
index represents a measure of the total control effort that must be 
expended to achieve the optimum DO profile. It is to be noted here 
that the values of the weighting factors are chosen prior to solving 
for u(t,x). 

In a similar paper Perlis and Cook [18] have used the same 
mathematical model as in [23] but modified to include two controls, one 
an artificial aeration control and the other an effluent control. The 


objective functional is correspondingly modified as, 


t x t x 
f £ ie 9 f ifs £ 9 
J= Ly J (VSP, ~ Vv) dxdt + Ly J if (U,) dxdt 
t x t x 
fo) fo) fe) fo) 
te Xe ; ce Xe . 
+ C, { if (VSP, = V,) dxdt + C, { ‘t (U,) dxdt 
c x t x 
fe) (o) ra) ro) 


(3-6) 
where VSE, is the specified BOD level 
VSP, is the specified DO level 
vy is the actual BOD level 
Vo is the actual DO level 
v, is the dumping control 
U, is the aeration control 


C. and C, are the weighting factors, 


Ly» bo» Oy 2 


pain _cioaigeagmpince ss 
oo Stan Dei! i alae al 
iothe Zoaaot cokes tennant i 
a eae - 
“to eeutsv off Jad a 


hei na , 
vi, vary 9 09°91 T 
aitlvfoe a2 eal Siidienie 


esse 63 ta ro ni 600 fon snare 
a a af se ye 
att ee, car mane rs eran 
| ieee - 


bial te pete nih a 1 ae a 


eae : ; tm ha ie 
> vie 
; ie Sie oe". ‘mr . d 
i), Sian? rae 
a ; hess a a : 


ro ae 


eid Ap -) ee ‘te 
Nie<s 


" 


30. 


Thus equation (3-6) is an expanded form of equation (3-5) to 
handle the problem of determining a management policy for the controlled 
dumping of effluents as well as artificial aeration. The values of the 
four weighting factors, Ly» Lo» Cy and Cy are chosen individually prior 
to proceeding with the minimization of the objective functional. 

3.3 Performance Criteria for Lumped Parameter Model 

In this section we shall formulate a general performance index 
for the lumped parameter state space model. This will reflect the various 
demands placed on a section of a river due to the addition of organic 
waste. In this task we shall use cost functionals proposed by earlier 
investigators for specific purposes. If a certain quantity of BOD is 
released from a waste treatment plant, there will be a cost associated 
with the processing of it. A typical relationship between the cost of 
treatment and the level of BOD expelled from a given plant is shown in 
figure 9(a). Such a relationship has been proposed by several investig- 
ators [6,10,15,25]. Another cost function relates the cost to the 
environment and the concentration of BOD. Obviously the former increases 
with the latter. A typical plot is shown in figure 9(b) [16]. The cost 
to the environment refers to the lost recreational facilities and other 
activities due to pollution. If we combine figures 9(a) and 9(b), the 
result is a function which represents the cost of processing waste to 
lower the level of BOD and the cost of excess BOD on the environment 
and the society (See figure 9(c)). Loucks [15] constructed a similar 
convex function for the target allocation of water versus the quantity 


of water desired by a user, while Bibbero [2] has used previous material 
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on cost-effective approach to establish the total societal cost for air 
pollution control. 
By fitting a second degree polynomial to this function in 


figure 9(c), we define the performance index for a section of a river, 


Ye 
= 2 
J = it (a, - a, u, + a, u,) dy (3-7) 
re) 
and a, = Fi 
| 
[2tr, +8, @,- 1! | 
ia | u | ee > (3-8) 
be m™m 
ae FE, + F, Cu, - 1) 
3 2 
u 
m 


where Fis Fy 


y is the distance along a section of a river 


and ui are defined in figure 8(c) 


Ye is the end of a section of a river 


u, is the BOD dumping control. 


ii 


It is useful to formulate some more performance indices, which 
can be combined with that given in (3-8), to form complex objective 
functionals. When using a section of a river for the dumping of waste 
the effects on the next section should be minimal. This can be taken 


care of using terminal costs given by, 
Se xy ea: (3-9) 
the level of BOD at the end of a section of a river, and 


J, = x9 (Y,) ) (3-10) 
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the difference between the saturation level of DO and the actual level 
of DO at end of a section of a river. 
3.4 Linear Combination Objective Functional 

Having formulated the three separate performance indices we 
shall next formulate a linear combination of the three indices. It will 
be shown in the next chapter that unlike earlier investigators, who used 
fixed values of the weighting factors, the optimum controls will be 
determined in terms of the weighting factors. 

It will also be shown in the next chapter that once the linear 
combination problem has been solved the results can be used for any 
objective functional expressed as a function of the given indices and 
the optimal controls can be determined. The conventional methods 
requires that the optimization problem be solved for each objective 
functional separately which drastically increases the amount of work 


to be performed. 
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CHAPTER (4) 
OPTIMIZATION OF WATER QUALITY 

4.1 Introduction 

Having derived an appropriate mathematical model in state space 
form for the river and having described the water quality management 
criteria by appropriate performance indices, the stage is now set for 
the optimization of water quality which is the subject of this chapter. 
In section 4.2, an introduction of the necessary control theory connected 
with the optimization of multi-cost system is provided. In section 4.3 a 
policy of water quality management is established for the controlled 
dumping of effluents (BOD) into a section of a river, subject to several 
different objective functionals. In section 4.4, the effects of artificial 
aeration is studied by introducing a vector control [uu] where u 


1 


represents the controlled dumping of effluents (BOD) and u, represents 


2 
the addition of dissolved oxygen (DO). The optimization is carried out 
with respect to the same objective functionals used earlier with uy alone 


as the control. 
4.2 Optimization Theory for Muiti-Cost Systems 

In this section the necessary theory related to the optimization 
of multi-cost systems will be summarized. For a more detailed discussion 
the reader is referred to the relevant literature [20]. 

Let the model of the system be represented by a set of n first 
order ordinary differential equations, 

x(y) = £ [x(y), u)] 


(4-1) 


x(0) = X 


«Echols ik» abbas A Salli Sissy. or Aaa 


Ballwaiaon 29 20}, hoddnidai8S 2 sadasgauke xaMlaup aosew,to xpatog 
\ePoV99 oF iosidila ,TevET = Ta) noly>06 e cAnk (G04) wtridulIe. 20 gabgewd | 


Letaktiame 46 nie} Ys oft don leser Gi .afsaoh yoni gviioetdo. aneawhaeb 
pi etede “Lat wl Wisden ada » stonhosant' yd hathugs oh cabanas 
ednonerye2 .u hire Coon} nana 20 satqoush belioszaen #3 sassesrgen 

tuo belyags ai fol jechebigo wot! (0a): nagyxo bev iosalh t© wotsthhm odd 
soins fa dtiw salizes been absnotdsrui evissatde aise of of tosquet dole 


.torsa09 od) ap 


336 


where f is a n-dimensional vector valued function continuous 


in x and u 


| 


is a n-dimensional vector representing the state of 


the system 


|e 


is a r-dimensional control vector of the system. 


The N performance indices are expressible in the form, 


a v ay % 
Jy = &elx (y,)] + i tix), uG)]dy, kel, wy N 
re) 
(4-2) 
or reformulated as 
2,(y) = &lx(y), uy) ] + (2g,[x(y)]/ax) - 
£[x(y), uly) ] 
(4-3) 
z, (0) rs g,[x(0)], k=l, ... , N 
Z-(¥-) = J, 
An objective functional is next defined as 
> = (2,5 ++ » 2y) (4-4) 


The control problem can be described by the following system 


of n + N = 1 equations, 
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[> 


. £, u) 


zxm % (xs u) + (98, (x)/9x]" + E(x, u) = G(x, w) 


6 = [86/8z,] - G(x, w) +... + [86/82] + G(x, u) 


and (4-5) 


2 (0) = e,[x(0)] 


=== 

Fam 
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= o[z, (0), eens Z,(0)] 


To obtain the necessary conditions for optimality for the free 
right end point problem, that is x(y¢) free, with the final distance, 


Yeo fixed, the Hamiltonian is defined as 
E 
H= P+ £@, uw) +P 6, @ w tee +P GG w 


Lap GG ut... +32 ac, »| (4-6) 


E 
nt+N+1 Oz 1 r) Zu 


where P is a n-dimensional vector. We get 


P = -0H/ox 
Bray 7 ~OH/dz,, kel, -.. 5 N (4-7) 
Paid Dey eo 


Using the fact that the right end point is free we have 
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P(y,) = 0 
Pek (Vg) a 0, k=l, eers N (4-8) 
Pent Y¢) adc e 


Now the expression for the Hamiltonian becomes 


i= Ee - £@, w) + @iy) 7 9¢/22,) ° GG, w) + eee 
+ @igy 7 96/940) * Goes a) (4-9) 


It has been shown [20] that the quantities Pak - 9$/dz,), 
k=1, ..., N, are constants and hence can be treated weighting factors 


Oye Using the second relationship of (4-8) the weighting factors become 


cy = ~39/3z, (4-10) 
Y Je 


The Hamiltonian can now be rewritten in a simpler form as 
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The procedure for the minimization of any objective functional, 
¢, consists of the following steps: 
Step (1) Consider the optimization problem with the same system of 
equations and conditions as in the given problem but with a 
linear combination objective functional, = 


k=1 
the optimum control, u*, as a function of the Cys. 


CyZye Determine 

Step (2) Search for the optimum values of the Cys which minimize the 
given ¢. This is accomplished by searching for the location 
where the gradient of ¢ is zero either with respect to the 
weighting factor ratios or the performance indices. In other 
words the search is made either in the cost scales space or 
the weighting factor space. 

Step (3 Determine the optimum control for the given > by substitution 
of the optimum weighting factor ratios from step (2) into the 
expression for the control obtained in step (1). 

4,3 Management of Water Quality by Controlled Dumping of Effluent (BOD) 
As our first attempt at optimization of water pollution, we 

shall consider the problem of determining the optimal dumping of effluents 

for the following three cost scales or performance indices. 
Let 
ae ey 

represent the level of BOD at the end of a section of a river. 
Ret 
J = X5(¥-) 


represent the difference between the saturation level of DO and the 
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actual level of DO at the end of a section of a river, and let 
Ve ; 
J, = if (a, - ayu, + aU, ) dy 
° 
represent the cost of processing waste to lower the level of BOD and 
the cost of excess BOD on the environment and to society, such as the 


aesthetic and recreational benefits. 


For the optimal control problem the system of equations is 


given by 
x (y) = -B x, (y) + 4, %) x, (0) = X15 | 
x (y) = C x, (y) - A x, (y) X, (0) = X59 
2,(y) = 0 2, (0) = X19 (4-13) 
zy (y) =O z, (0) = Xo 
za(y) = a, - au) + a,luj@)]? — 2,(0) = 0 


We have to satisfy the following conditions: 
(i) final conditions, x(V_)> are not fixed 
(ii) the dumping control, u(y) is unconstrained 
(iii) the distance, Yeo is fixed 


We shall consider the following objective functionals: 


(4) (2,5 2s 23) = 23 > (c,/e,) 2) 
ta (c,/c,) Zo (4-14) 


(41) $(2), 25. 23) = 4324 * 29 + 25 (4-15) 
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(411) $(2,5 255 2,) = a, (2, + 2,°)/2 +2 (4-16) 
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Solution of the Linear Combination Problem 
The first step, namely, the solution of the linear combination 
problem is common for above three objective functionals. Using the 
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near combination objective functional, kal 929 nea. 


Hamiltonian becomes 


aS Z 
H c, (a, Agu, + 2g0, )- Bx,P) + u4P, + 
Cx,P, - Ax,P5 (4-17) 
And 
oH 
By Suaciw eR iene Pas P; Ve) He, (4-18) 
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alee os Ox, AP, ; Py (Y¢) Co (4-19) 


After solving the costate equations, the expression for the 


costate variables are 
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With no constraints on the control function the optimum 


dumping control takes on the following form, 
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fac CCA wD 
1 1 G | aD | 
where Ke == | |_| - eel | 
ss 2a, \e3 (B-A) | a 
2 2a, | (B=A) Dy) 


Substituting the optimum dumping control back into the system 


of equations (4-13), the BOD profile is given by 


a —-By B(y-y,.) _-Bly-y-) 
“ By (peepee ile eee a ee): 
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(B+A) (4-23) 


= Ky 


The profile for the saturation level minus the actual level of 


dissolved oxygen is described by the expression 
—Ay —Ay Pale 
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Using equations (4-23) and (4-24) we can obtain the expressions 


for the performance indices as, 
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(4-26) 


-2Ay 
2 (le f) 
+ Ky TA (4-27) 


At this point, in order to point out the optimization 


procedure effectively we shall use the numerical values used by Perlis 


and Cook [18]. Let 


K 


4 7 0-66/day 


Ka K. = 0.16/day (4-28) 
V = 1.00 miles/day 


In addition to the above values, the distance under consideration 


is chosen equal to two and one-half miles. Using these values, the three 


performance indices become 


f 


1 
Z, = 0.6703 x19 + a; nee a, - 0.8605 


1 (4-29) 
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as defined in section 3.3. 


In reality, the cost of processing waste, Fie approaches 


infinity as the BOD concentration is lowered to zero. This allows us 


to rewrite the above as 
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Minimization of the First Objective Functional 


Substituting the expression for z z, and z, from (4-32) thru 
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To determine the optimum values of the weighting factor ratios, 
c,/c, and Co/Cas which will minimize the objective functional, we take 


the gradient of 4 and set it equal to zero to get 
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Checking to see if the weighting factor ratios (4-36) and 


(4-37) minimize or maximize the objective functional, we used the second 


derivative test. We have 
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The objective functional is minimized for this particular set of weighting 
factors. 

Substituting the optimum values of the weighting factor ratios 
into the expressions for the performance indices and the objective 
functionals, we get their minimum values in terms of the initial conditions, 


x and x and the BOD level criterion, ui for the dumping of organic 
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After carrying out the appropriate calculations, we are able 
to establish the following relationship: 
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and x, for realistic values of 


To determine the shape of Uys Xy 9 


X10? *o0 and ui» we use the same values as Perlis and Cook [18], namely 
X19 = 30 mg/2, Xo = 3.06 mg/2, and wt 23 mg/2, to evaluate the system 
of equations. The results for the dumping control, the level of BOD and 
the saturation level of DO minus the actual level of DO are plotted in 
figure 10. The particular values taken on by the weighting factor ratios, 
the optimum control, the performance indices and the minimum objective 


functional are furnished below. 
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Minimization of the Second Objective Functional 


In a similar manner, we proceed to determine the optimum 


dumping control for the second objective functional, ¢ = a, z, + z,+2z 


Spl 3 3; 
given in (4-15). As outlined in a previous section, there is no need to 
solve this problem from the beginning. Since the objective functional 


fs made up of the same cost scales z and z,, the solution to the 


Tree 3 
linear combination functional earlier can be used for the second objective 
functional. Substituting the performance indices (4-32), (4-33) and 


(4-34) into the new relationship of 4, we can write, 
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As before the gradient of » is determined and set equal to 


zero and the ensuing equations are solved to give the following optimum 


values. 
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Using the second derivative test, we can establish whether 
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which fulfills the conditions for the objective functional to be a 
minimum for this set of the weighting factor ratios. We can evaluate 
the relationships of the performance indices and the objective functional 


to determine what their minimum values are for this case. 


(May eee, og? 1000+ salah = (ohgo? 
C= gt i SON ge PELIT ~ ye HO > fyo\g2) 


Yoilterhe inal: OW .998' or sbi tea s83 gatel 
ep) ‘nso Sviisstdo ots setnaat (200) bra 


ee oY 


ent me 


rrigainge 


503 


CF 0.5046 X19 7 0.1368 Xoo t 1.6957 us (4-46) 
Z» = 0.1309 X19 + 0.1758 Xoo tb 0.2165 Us (4-47) 
z, = a.,{0.01645 x a 0.01092 x ‘ + 0.08281 ‘ 

3 gies 10 ; 20 : un 


+ 0.02610 X41 Xoo + 0.07353 X10 us 
+ 0.05715 x59 u,] (4-48) 
Mid i eye Sey Sy aS 
310. 10 : 20 . i" 
+ 0.09687 X10 Xo + 0.4048 X19 un, 
+ 0.3256 x9 u,] (4-49) 


Using the same values for x and uu the optimum profiles 
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are plotted in figure 11. The values of the optimum weighting factor 


ratios, the optimum control, the performance indices, and the minimum 


objective functional are 
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Min 96 = 623.2671 a, 


From these results we can conclude that the optimum weighting 
factor ratio, c,/c., is approximately equal to the performance index, 
Zo» multiplied by ay, while the other optimum ratio, Co/Cas is 
approximately equal to A321. 
Minimization of the Third Objective Functional 


We have 
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Again we notice a direct relationship between the optimum 
weighting factor ratios and the two performance indices to state that 


the ratio c,/c, approximately equals a,Z4 and the ratio Co/e, approximately 


equals AxZ5- 


4.4 Management of Water Quality Controlled Dumping (BOD) and Artificial 
Aeration (DO) 
In this section, let us consider adding an aeration control to 
the same section of the river to see what effect it has on the profiles 
BOD and M0 gar - DO). Now the differential equations, (2-7) and (2-8), 


for the system are written as 
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where u,' is the dumping control 


Uy is the aeration control 
As in the previous case, the state variables are given by 
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The performance indices, Jy and Jos will remain the same but 


J, will be changed to accommodate the second control, Uy. 
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We consider first the linear combination of these indices as objective 


dy (4-53) 


functional. 
Solution to the Linear Combination Problem 


The augmented system of equations is given by 
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The costate equations and the solutions to them are the same as before. 
There are no constraints on either control, hence the optimum dumping 


control is 
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and the optimum aeration control is 
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Substituting the optimum controls back into the system of 


equations (4-54), the optimum BOD profile is 
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Using the same values as given by (4-28), the three performance 


indices are 
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Having solved the linear combination problem, we now optimize 
the system with respect to the three objective functionals as given by 
the equations (4-63) thru (4-65). 

Minimization of the First Objective Functional 


We have 
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The optimum weighting factors are determined in the same manner 
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To illustrate the optimum profiles (see figure 13) for this 
objective function we choose the ratio of a,/b to be equal to unity. 
The other parameters remain the same as in previous examples. The 
optimum weighting factors, the resulting expression for the performance 


indices and the objective functional are 
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2, = 0.4469 x, + 0.0000 x). + 1.3738 u, (4-71) 
z, = 0.1021 x), + 0.1280 x54 + 0.1784 u, (4272) 
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2 2 2 
Zz = [0.02977 X19 + 0.005558 Xoq + 0,2742 uF 


+ 0.004113 x + 00,1787 x UT 


10 “20 10 


+ 0,.0008919 Xo uw a., (4-73) 
Minte dom OROA0) date 01667 me ce a C2 
J 10 ¥ 20 F m 
+ 0.01234 X19 *20 + 0.5360 X19 Un + 0.002684 XQ ui a, 
(4-74) 
For Xj = 30 mg/2, Xo = 3.06 mg/% and Use i 23 mg/2 
(c,/c) = 25.8070 a, 
(co/e,) = 2.7812 ay 
u, = 23. - 13.3485 exp [.16(y-2.5)] + 0.4450 exp [.66(y-2.5)] 
z= 45.0030 
Zz, = 7.5569 
2, = 295.6123 a. 


Min > = -886.7966 ay 


We notice that two of the performance indices are related to 


the weighting factor ratios by 
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2a Piao 0517214 c,/c, 


a (4-75) 


= | 
: ; Poets “Pood Cy/cy 
To compare the optimum solutions of the state variables for the 
system with the aeration control and without this control, the optimum 
trajectories of the state vector, [x, ihr are shown for the section of 
the river under consideration in the figure 14 on the XX plane. From 
the plot for the state plane, we find the initial and final conditions 
of the trajectories are the same for both problems when the optimum controls 
are applied to the system. 
Minimization of the Second Objective Functional 


The second objective functional, $= a24°Z5+Z3, for the case of 


the aeration control takes on minimum value when 


(c,/c,) = [-0.09677 X10 + 0,.2375" x, 07 0.04686 ui ay (4-76) 


20 
(co/c) = [0.6818 X19 7 0.1848 X59 Te2e2915 ui a, (4-77) 
For X19 = 30 mg/2, X59 = 3.06 mg/% and dee 23 mg/2, we plot 


the optimum dumping control, the optimum aeration control and the resulting 
State variable profiles in figure 15. In addition the weighting factor 
ratios, optimum controls, performance indices and minimum objective 


functional becomes 


(c,/c,) -14.7932 a, 


72.5909 a 


(co/c,) 3 


u, = 23. - 4.2179 exp [.16(y-2.5)] + 11.6145 exp [.66(y-2.5)] 


and vol esidaiany oteis ont 
momiago sii ,Lors0os atda : STIMOD aM 
te woliose ads ‘rot nworla os Jt i . 2398 eT 
novi wang ox ,x sda 0) bE mavig?2 ods at #ottetsbtes 
anolitbmos Isat? bos tigen oft batt ow «iat 2820 sid x02 solg ods 
alozinos awmitqo eft notw onal dong, Ajod s0% sane ‘sini eaiaesei tere ons x» 7 
ee sien betiqqs sxe _ 


40 9882 913 Tok; Shes tee -tsaotioayt 1 rte aan . —_ 
'! pedw suey muimtabm wo tesley Lo1Jno0 aokdmson sit 
7 

(BV-8) -8 [ou 38R2.0 = gx 2TESO + Gx TNOEDIO-). = (9\z9) 


(vt=a) e® [uv S10S.S + ocx BRBL.O = or 8183.0) - (p2\go) 
tolg ow -t\ga CS = bas tge B08 = go: 8 \ ger OF gx 20 (= 
gottioser edi baa. Loriaos dottnaun taut 2¢0 sit Toxsige. gakiauh matzgo sid 
403082 adtsdgbew ods nakatbbs ‘at «2S ozvgtt Poon hreabe a - 
avissstdo mmintm bre am cass mothoqo ,eolts: 
| swans Lembo t 
ae : 
- 


[(BaS-¥)02. i ehnp.L nieieon oo as ase 77 


| 


te" 
‘5 
> “ 
, Ss by a : : 


7 a ‘Yy 7 aa, a fo, Pare oe — "Aa, 


63: 


f without 
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FIGURE 14: State Plane for ore, ~ (c,/c,)2, - (c4/c,)2, 
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uy, = 36.2955 exp [.66(y-2.5)] 
Z, = 72.5885 

Zo = 14.7975 

Z, = 986.1277 a, 


Min 9» = -87.9966 a. 


The relationship between the optimum 
and the indices, 24 and Zo» remain the same as 
the aeration control, but the actual numerical 
In figure 16, the state plane trajectories are 


and without the aeration control. It is to be 


diverge from one another. 


5 


ratios, c,/c., and Coleg, 
the previous problem without 
values are not the same. 
plotted for the system with 


noted that the trajectories 


Minimization of the Third Objective Functional 


For the third objective functional, pra, (2, +2,-)/2 + Zz 


get the optimum weighting factor ratios as, 


(c,/c,)= [0.3555 x, - 0.007901 x9 + 1.1013 uJ a, (4-78) 


20 


(c,/c4)= [0.08372 x) + 0.1396 Xy + 


0 0 


for the a/b equal to unity. 


0.1093 ui! a, (4-79) 


Again we plot the optimum controls and resulting state variables 


for X19 = 30 mg/2, Xoq = 3.06 mg/% and ues 23 


evaluate the following, 


(c,/c,) = 35.9709 a. 


mg/2 (see figure 17), plus 
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(cy/c4) = 5.4519 a. 


u, = 23. — 18.8578 exp [.16(y-2.5)] + 0.8723 exp [.66(y-2.5) ] 


uy = 2.7260 exp [.66(y-2.5)] 
Zz = B5e9157 

Zy = 5.4679 

Zz, = 583.0876 a, 


Min 9 = 1245.1621 a. 


As in the case of the system with the dumping control only, 
we have the optimum weighting factor ratios equal to the performance 


indices multiplied by a Again plotting the trajectories for the state 


3° 
plane in figure 18, we are able to conclude that the trajectories are 
similar. In other words one can be considered to be a reduced mirror 
image of the other and the final conditions of the state variable, X19 
are nearly equal in magnitude for both cases. 
4.5 Sensitivity Analysis 

In this section we will present a discussion of the sensitivity 
of the objective functionals used to changes in the weighting factor 
ratios. Such an analysis yields information about the behaviour of the 
objective functionals in the neighborhood of their minimum for small 


deviations in the value of the optimum weighting factor ratios. A 


quantity called the objective functional sensitivity vector is defined by 
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9 { 

o(e4/ cep | 
22 a "3 | , nal,2 (4-80) 
=| C@erry C) 

3(c,/c,) | 


| c,/c, = (c,/c,)*+8 


where cafe,/c, is the weighting factor vector 
c,/c, 


(c,/c,)* is the value of the optimum weighting factor 
ratio 
A is the deviation from the value of the optimum 


weighting factors 


This definition is similar to the performance index sensitivity 
vector defined by Sage [21]. It is the gradient with respect to the 
weighting factors evaluated at c = c*¥ + A. Now, the sensitivity of a 


objective functional can be expressed as 


ah 
ck} 9 
25 PY Bae Shy d9(c,/c.) et 9(c,/c,) A 
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c=cr*+A af (@ ° 
Pee hfe —=— +A phe 
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Rae 


(4-81) 


The values off the sensitivity and sensitivity vector are 


listed in table 3 for each objective functional. 
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From the table, it is seen that the first objective functional 
is most sensitive to the deyiation of the weighting factors from their 
optimum value while the second objective function is the least sensitive 
in both cases without the aeration control and with the aeration control. 
In general for all three objective functionals, the sensitivity vectors 
are most sensitive to the deviation in the weighting factor, cy/Ca, 
than the weighting factor, Co/ca. 

4.6 Interaction of Performance Indices 


Next, we shall relate the performance indices, z and z 


120%) 3? 


to each other to investigate how they vary with regards to each other 
for both cases: without and with aeration. Turning our attention to the 


equation for the performance indices, (4-32) thru (4-34) and (4-63) thru 


(4-65), we notice that the three eqyations are functions of two variables, 


/ 


c,/c, and Cy/ca. This allows us to combine the equations into single 


equation while eliminating the weighting factors from it. Using the 


same parameter as before, the equation for system I becomes 


2 


2 
Z,/a. = 2.5660 Z, 32.4373 z ATS 2 69220 Zo 


Te? 


qr Pal Rae) Zz 7 814.7469. .2,. + 3900.2377 (4-82) 


iz 


and the equation for system II becomes 


vd 


2 
z,/a, = 0.6014 2, = 0.3319 21 2» + 1.3560 Zo 


~ 77.4346 Zz. - 8.3379 Zo + 2660.5918 (4-83) 
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Both equations are second degree in Zy9 Bo and Za To determine 
the classification of these quadric surfaces, a rotation and translation 
or coordinates in the three-dimensional space will be performed. The 
rotation of coordinates in the plane, 25 = ucos0—-vsiné and Zo = usinO+vcos6, 
is used to dispose of the Z1 25 term in each equation. Rotating the 


relationships by 7.005° and 11.872° for system I and system II, respectively, 


the equations of the surface becomes 


2 


z,/a, = 0.5725 uv + 134.5154 fe AST GaN 


— 811.2579 v + 3900. 2377 (4-84) 


z,/a, = 0.5665 ae rt) Sameer eckia 


Tome OO8V t+ a2000.5918 (4-85) 


By introducing the translation of coordinates, x = u-68.3496, 
y = v-3.0155 and z = Z,/a,-2.4767, we find, for the equation of the 
surface of system I, 
pee ie + inks an mz (4-86) 
(1.3216)* (0.08622) 
which we recognize as an elliptic paraboloid. 
The translation of coordinates, x = u-68.3968, y = vt2.7934, and 


z= Z4/a,+0.4178, shows that the surface of the performance indices for 


system II is also an elliptic paraboloid with the equation 


2 Z 
eee et ao Jat ES = Z (4-87) 


(1.3286)7 (0.8479) 


‘ - . 
wee’ 


ee ee 


: eu ; -_ 
a i - * se dai ° Su actea 7 e*\y* 
(48=0) | “VES .OORE + -— Of2o.E8e = 9 
Vig a 
. <6 age 7 
(28=s) a BLO2.089S + WOOT Te | eT 


wv d6e8,t0 = *y covet + 8, eada.o = gales 


JCHE.BO-u = x ,aeIAatbIOd Yo motrmLenexs odd yrlovborsat ye | NT 
sd to nolssups: sd) 702, S829 SH ,\aTO.S-pe\,x = 2 ns 2et0.¢29'= y 7 


74. 


The two surfaces are shown in figure 19 for the new coordinate 


Systems. The mapping equations are 


System I: x = z, cos 7.005° + z, sin 7.005° - 68.3496 
y = -z, sin 7.005° + Z cos 7.005° - 3.0155 
z= Z,/a. - 2.4767 

System II: x = Z, cos 1198 72- ot Zo sin 11.872° - 68.3968 


y= 2, sin 11,8722 tez" cos) 11.8720 42257030 


Z 
z= Z,/a., + 0.4178 


with the inverse transformations given by 


System I: z, = (x + 68.3496) cos 7e00o: te (y + 3.0155) sin 7.005° 


1 
2) = (x + 68.3496) sin 7.005° + (y + 3.0155) cos 7.005° 
Z, = (z + 2.4767) a. 

System II: z, = (x + 68.3938) cos M872 (y= T7984) cin 11.872 — 
Zz = (x + 68.3968) sin 11.872° +207 — 2.7034) fos 1f.372- 
Zz, = (z - 0.4178) a, 


From figure 19, the surface for the system with aeration control 
has a very much larger minor axis than the surface for the system without 
aeration control while the major axis is approximately the same for both 
cases. All three optimum solutions of the objective functionals of 
figure 19 a) appear quite close to the major axis in the back quadrants, 
and only two optimum solutions of figure 19 b) appear in the same region 
with the third solution corresponding to the functional, 9 = a, Zz ° 25 + Ze, 


near the minor axis. For the system without aeration control the solution 
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of the functional, 9 = a, Zz. Zo + Z3 is the closest to the local 
minimum of the performance index surface while the functional, 

> = ay ae + z5°)/2 + Zs is the furthest away. The order is changed 
for the system with the aeration control, since the optimum solution of 
¢ = 24 - (c,/c,) Zz. - (c,/c) Za is now the closest and 4 = a, 2, * 2 + 2, 
is the furthest away. The significance of these observations is to 

provide a means of determining the best objective functional by considering 


the location in relation to the local minimum of the optimum solution 
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CHAPTER (5) 
CONCLUSIONS 

5.1 Summary of Results 

In the previous chapter, we determined the optimal water quality 
Management policy for a section of a river using multi-cost contnol 
techniques. Two cases were considered: first, with one control, namely 
the dumping of effluents (BOD) only, and second, with two controls, the 
dumping of effluents (BOD) and the addition of artificial aeration (DO). 


The results obtained are summarized below: 


I. System with One Control 


(a) System Model 


x, &) egies eee | 1 
| ie 
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x, (y) 
where A, B and C are defined in (2-11) 


(b) Performance Indices 


(i) J, = x, (Y¢) 


(it) J, = x, (¥¢) 
7 2 
(1ii) J, = J (a, “a, Uy + a, Uy ) dy 
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(c) Objective Functionals Used 


(i) o# oe (c,/c,) 2,7 (cy/c,) Z5 


(ii) 9 = az 24°25 +. 2. 
@,° + z5°) 
(iii) 9 = a ape I Z. 


(d) Optimum Controls for Each Objective Functional 


(i) u, = 23.0 - 39.9exp[.16(y-2.5)] + 43.4exp[.66(y-2.5) ] 
(ii) u, = 23.0 - 13.3exp[.16(y-2.5)] + 8.6exp[.66(y-2.5) ] 
(iii) u, = 23.0 - 19.lexp[.16(y-2.5)] + 1.2exp[.66(y-2.5)] 


IME System with Two Controls 
(a) System Model 
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(b) Performance Indices 
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(c) Objective Functionals Used 


@) o= Zz - (c,/c4) Z4 (c,/c,) Z. 
(ii) 6 = a, 2) Zo + Z., 
ye 2 
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(d) Optimum Controls for Each Objective Functional 


iby, oe 23.0-13.3 exp[.16(y-2.5)]+0.45 exp[ .66(y-2.5)] 
u, = 1.39 exp[.66(y-2.5)] 

(ii) u, = 23.0-4.2 exp[.16(y-2).5).]+11. 6 exp[..66(y=2.5),] 
Uy = 36.3 exp[.66(y-2.5) ] 

(iii)u, = 23.0-18.9 exp[.16(y-2.5) ]+0.87 exp[.66(y-2.5) ] 
uy = 2.73 exp[.66(y-2.5)] 


5.2 Discussion of Results 
A study of the results shows that the optimum dumping controls for 
the first two objective functionals dump more organic waste in the second 
half of the river section. The opposite is true for the dumping controls 
of the third functional. This is due to the increased weighting placed 
on the final conditions in the third objective functional. The aeration 
controls are the same for all the functionals and differ only in magnitude. 
Now the dumping and aeration control problems of Chapter 4 can 
be compared to work done by Perlis and Cook [18]. The optimum controls 
which Perlis and Cook calculated according to their objective functional 


(see equation (3-6)) are quite different from the ones determined here 
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for the three objective functionals given by (4-14), (4-15) and (4-16). 
Their dumping control has an initial and a final value of zero with the 
majority of waste released in the latter half of the river section. The 
BOD profile decreases in magnitude along the river section which is 
opposite to the increasing BOD profiles determined in the previous 
chapter. Their aeraticn control has a non-zero initial value and a zero 
final value with the majority of oxygen released in the first half of 
the river section. The DO profile obtained by them are similar to the 
results obtained in Chapter 4. Both have the same concave upward 
appearance. 

The optimization theory for multi-cost system in Chapter 4 has 
given an overall optimum control policy for the with and without aeration 
cases. This method of arriving at an optimum eonerel policy for the 
management of the water quality has proven to be more flexible in its 
approach then previous papers on the subject. The results have contributed 
to better understanding of the behaviour of the system for different 
objective functionals through the investigation carried out in the 
previous chapter. 

5.3 Suggestions for Further Research 

The discussion in this thesis is not intended to set down a 
comprehensive Water Quality Management Policy; what has been attempted 
is the outlining of a procedure to be followed to obtain one such policy. 
We have shown some of the possible results that might be arrived at by 
using multi-cost control techniques. The following are some areas for 


further work. 
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In Chapter 2, a bilinear model was introduced for comparison 
with the lumped parameter model, which showed close agreement 
over the limited range investigated. A rigorous study 

could be conducted in this area, followed by the application 
of multi-cost system theory to determine the necessary 
optimum controls. 

The work carried out in Chapter 4 can be extended further 


by the introduction of some addition objective functionals 


and the determination of the resulting optimum weighting 
factor ratios. 

In section 4.6 of. Chapter 4, it was shown that thee n” 
performance indices can be related to one another by 


"n-1" weighting factor ratios to give a 


eliminating the 
surface. Each point of this performance index surface 
corresponds to a solution for an objective functional. 
When the surface has a local minimum, there would be an 
optimum cin emanate Using this approach, the multi-cost 


theory could be extended without knowing or defining 


optimum objective functional. 
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APPENDIX I 


Solution of Homogeneous Equations by the Analog Computer 


System ig 
“ Kr 
B=- V B 
‘ Kd Ka Ka 
=-—B-~D+—~—D 
2 V } V V 


LC. 
D 
sl ey 
\ 
Parameters Used 
K = 0.00 > 0.16 1/day Dg = 9,00 meg/2 D(O) = 6.00 mz/% 
2 

Ka ="0, 1671/ day Vi =) 00) mee B(O) = 15.00 mg/2 
Ke = 0.66 1/day 
Results 


See figure (4). 
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System Ze 


D=- Boe Dt ay) 


Computer Diagram for the Solution of System 2 
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Multiplier B 
AK 
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<| 


Parameters Used 


K = 0.00 > 0.05 2/day*mg Daa 9.00 mg/2 D(O) = 6.00 mg/2 
Ky = 0-06) 1 day V = 1.00 mile/day BCG), = 15.00 me/2 
Maeroco L/day 

a 
Results 


See figure (5). 
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APPENDIX IT 


Comparison of Lumped Parameter and Bilinear Models 


System ies 
e K.. 
Berea 
4 Ka La K 
Ve ayaa & 21 ta 2s 
System 2: 
‘ K 
Be wie 7? 2 
C Ka Ko K 
Vee haya ot eee 
Parameters Used 
Ke = variable K = LAK D, = 9.00 me/2 D(O) = 6.00 mg/2 


i} 


15.00 mg/2 


K, = variable K.= variable V 1.00 mile/day B(O) 


ther BO.a =~ (Ou = WAgm O08 = Sab eee 
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